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Abstract 

The relativistic mean field model, the Zimanyi - Moszkowski (ZM) Lagrangian describes 
nuclear matter and stable finite nuclei even in the non- relativistic limit. It fails, however, 
to predict the correct non-relativistic spin-orbit (SO) coupling. In this paper we improve 
on this matter by an additional tensor coupling analogous to the anomalous gyromagnetic 
ratio. It can be adjusted to describe the SO-term without changing the mean field solution 
of the ZM-Lagrangian for nuclear matter. 
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Introduction 



The advantages of describing nuclear physics relativistically are to find connections between 
different parts of the non-relativistic mean field nuclear potential, such as the linear en- 
ergy dependence emerging from kinematic features of an effective Dirac equation for single 
nucleons or a mild mean potential being balanced between strong attraction and repulsion. 

The most known such model, the Walecka model[|l], 0] was first intended to be an elementary 
description, but it was soon discovered that only an effective Lagrangian approach is fruitful. 
Renormalization problems and divergences of higher meson loop corrections with a strong 
effective coupling revealed that by dropping the requirement of an elementary approach a 
higher degree of freedom in the choice of the effective Lagrangian can be exploited. 

The Walecka model describes the non-relativistic mean field potential, (V — S), and spin- 
orbit coupling J^(U + S), both derived from the same effective Dirac equation. Unfortu- 
nately it fails to describe compressibility and a finite nucleus, which remains bound in a 
quantum molecular dynamics (QMD) approach. Simulation of dynamics including heavy 
ions reactions at relativistic energies, however, desperately needs a Lagrangian, which is 
relativistic and has a stable nucleus. Such an effective Lagrangian was proposed by Zimanyi 
and Moszkowski || recently. It works with less repulsion, V, with an attractive mean field, 
S which is bounded from above by the free nucleon mass, M and describes stable finite 
nuclei |||-[§|. Up to now only the too weak SO-potential was a manque of this model 0. 

Recent investigations did show that this Lagrangian has many advantageous properties in 
describing the bulk properties of the nuclear matter |T0[|-||12||. The optical potential derived 



from the ZM Lagrangian also agrees with experimental findings fl3H. Finite nuclei described 



within the framework of this model are used in QMD calculations! 14] as well. 



In this letter we show that an additional tensor term influences the spin - orbit part of the 
non-relativistic potential without spoiling the benevolent property of the ZM Lagrangian: 
the mild repulsive mean field. This reduced repulsion is necessary to describe stable finite 
nuclei in the quantum molecular dynamical approach to heavy ion collisions. Throughout 
this paper we use Minkowskian metric with the signature (+, — , — , — ) and the convention 



1 



h = c = i. 



Improved effective Lagrangian 

The Zimanyi-Moszkowski (ZM) Lagrangian, Czm, is formulated in terms of nucleon fields 
■if), a scalar o and vector uo^ mean fields. The corresponding scalar and vector coupling 
constants are g s and g v , respectively. Now we add a tensor term to this Lagrangian which 
couples the antisymmetric spin current density ipa^tp to the field strength of the vector 
field F^ v 

£zmb = -^M^+^{[id^-g v u^)-^ + ^g v a^F^ij 

-\f^ + \m 2 u^ + i«V • &>a - U(a). (f ) 

with 

M 

M* = —— (2) 

l + g s a/M y) 

being the effective nucleon mass. M and m are the bare nucleon and vector meson masses, 
respectively, and U (er) includes the \m 2 a a 2 mass term and an eventual self interaction. 

In order to deal with static mean fields over a finite spatial extension we assume that only 
a and = (cuo, 0) differs from zero and all time derivative vanishes. In this description of 
static nuclei uq and a may have gradients, however. The field strength F^ v has only electric 
components F^ = — VfcC^o- The corresponding coefficient Dirac- matrix in the tensor term is 
a k = —ictk- The equations of motion for the meson fields can be derived from the improved 
Lagrangian easily: 

- Ha - U'(a) + ^ 2 gM + A-^Lg^a^ = 0, 

+ = 0. (3) 

These equations lead to the same mean field equations as without the tensor term for all 
unpolarized media with (ipa^ip) = 0. Using this the homogeneous static mean field solution 
arises as in the case of the Zimanyi - Moszkowski Lagrangian with U(a) = \m 1 (y o 2 : 

M * m\ 
171(7 a = Jp 9s ^'> 

m 2 uj = g v (^H)- (4) 



The Dirac equation derived from the above Lagrangian on the other hand has a new 
contribution 



-f(id„ - g^) + ^g^F^ - M, U = 0. 



(5) 



In order to investigate the non-relativistic limit we consider the Hamiltonian, expressing 
Hip = id t ip from the above equation. Multiplication with 7 from the left and neglecting 
vanishing mean field and time derivative components gives the relativistic Hamiltonian 



H = -iaV + g v u - -^rtfVg v u + M*/3, 



(6) 



where we used the Dirac representation of matrices 
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(7) 



For the sake of a more compact notation we introduce S = g s a / (1 + g s a /M) and V = g v cu . 
This notation automatically satisfies the relation M* = M—S. The relativistic Hamiltonian 
equals 



H = (3M* + ap + V - i t(W) 



(8) 



where the momentum operator p = —iV is introduced. The Dirac matrices and spinors 
in 2 x 2 notation, if) = (tp, x), lead to coupled equations for the upper (tp) and lower (x) 
components. The non-relativistic Hamiltonian is the effective Hamiltonian acting on the 
upper component of the Dirac spinor. We get the following eigenvalue equations 



(V + M*-E)v + t(p-i^ f Vvjx = 
a(p + i^ f VV\ l p + (V-M*-E) X = 0, 



(9) 



with E being the relativistic energy eigenvalue. The non-relativistic energy eigenvalue is 
therefore e = E — M. Introducing further simplifying notations, 

A 



p ± = p zp i 



2M 



(vn 



V± = V± M* — E, 



(10) 



it reduces to 



V+p + vp+X = 0, 
<?p_ + V-x = 0- 
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Here V± can be expressed by using the non-relativistic energy eigenvalue as V + = V — S — e 
and V- = V — S — e — 2M*. Extracting now x from the second equation and substituting 
back into the first gives 

(v + - ap + ^-ap^J = 0. (12) 

A non-relativistic Hamiltonian can be obtained by expressing the energy eigenvalue e from 
this equation. We get 

hip = eip = — S — dp + —ap^j <p. (13) 

The non-relativistic expansion of this operator together with a gradient expansion for the 
mean field factors V and S leads to the non-relativistic nuclear potential. The spin-orbit 
coupling term arises as a sub-leading correction in an 0(1/ M*) expansion. Traditionally it 
is due to the fact that the p gradient operator does not commute with 1/V-. In the present 
case we have an additional contribution to this term stemming from the tensor coupling 
due to p + x p_ ^ 0. We proceed by calculating 

P + ± - ±P + - l9±. (14) 

After noting that = V(V + S) both terms on the right hand side can be expanded in 
inverse powers of the effective mass. We get 

P + Y_ = ~ + ([e ~ (V ~ S)] p + iV [(1 + X)V + S]) + ... (15) 

Let us denote this result by — 7T+/2M*. Using this notation the non-relativistic Hamiltonian 
is 

h = V-S + ^a.n + a.p_ + ... (16) 

with 

7f + = p- -L ( [e - (V - S)] p + iV [(1 + A)V + S]) (17) 

Now we utilize the general algebraic identity, 

a • 7r + a • p- = 7r + ■ p- + io (n + x p_) . (18) 
The individual terms are again expanded 



and 



n + xp_ = ^V[(l + 2\)V + S] xp + x W + ... (20) 



The non-relativistic energy eigenvalue of fll6|) in an 1/M* is given by 



' h ° = v - s + &, + - < 21 > 

Substituting this leading order result into the coefficient of p 2 in the next order term we 
make a higher order error only. This arises as a relativistic correction to the kinetic energy 

1 e -(y-ffu 2 p 4 

fel = -2A£ 2M. P =~8M! + --- (22) 
Two further interesting term arise from the scalar product 7? + • p_: 

' h ^-wi^ v+s) ^ + ^yt)^ v - < 23) 

Factorizing the 1 / 4Mjf coefficient we get 

^ = -jk^ {v + syt - XM ^(if, tv ))- (24) 

This will be the analogue of the so-called Darwin term. After symmetrization and integra- 
tion by part we arrive at 

62 = 8W 2 ^ [(1 + 2X)V + S] + 4^P {1 + X)fS ' f{V + S) (25) 
The familiar non relativistic Darwin term stems from a Coulomb potential (V = Vc, 5 = 
and M* = M). In this case V 2 Vc oc 8(f) and it can be interpreted as the spin-orbit energy 
shift to zero angular momentum states. 

Finally the term stemming from the vector product between n + and p- leads to 

k = (V((l + 2A)K + S) x p) + (vA x w). (26) 

Acting by the gradient V on the 1/M* coefficient a result higher order in 1/M* arises, 
which can be consequently neglected. Finally we arrive at 

h = ^ (V((l + 2\)V + S)xp)+... (27) 

This result - and its derivation - is analogous to the treatment of an anomalous gyromag- 



netic ratio (A = Ag/A,g v = e), |TH] 



For central mean fields V = V(r ), S = S(r) the first term of this expression can be rewritten 
by using the angular momentum operator, L = r x p: 



Summarizing the main effect of introducing the direct tensor coupling term in the relativis- 
ts Lagrangian is to enhance the role of the vector mean field V to (1 + 2X)V in leading 
order spin-orbit corrections to the non-relativistic single nucleon energy. From this result 
we conclude that the effective strength of the spin-orbit coupling is increased by the new 
term in the Lagrangian, Vso — (1 + 2A)V + S. This is the main result of this paper. 

The equations of motion for the mean fields do include an additional term due to the tensor 
coupling. This has no effect, however, in the mean field approximation to unpolarized 
media. Therefore estimating the spin-orbit potential we can use the values of V and S 
obtained self consistently for the ZM Lagrangian (V = 82 MeV, 5* = 140 MeV). Comparing 
the strength of the respective terms with that of the Walecka model (V + S = 785 MeV) 
we obtain A = 3.5. One can improve on this estimate by applying the presently described 
Lagrangian, £zmb, for finite nuclei. 

It is interesting to note that P.G. Reinhard and his collaborators also considered a tensorial 
coupling in connection with the original Walecka Lagrangian. They discarded, however, 
this possibility in that case, because they found its contribution negligible^]. 
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